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ON '.T5 •■'•‘“-I BLUING IN A VISCOUS 
DIFFUSIVE STRATIFIED FLOW 

C. G. Koop and L. G. Redekopp 

Department of Aerospuuc Engineering 
University of Southern California 
Los Angeles, California 90007 

The effect of diffusion of specie upon the flow 
about a transverse flat plate moving horizontally 
In a viscous stratified medium is considered. 

Asymptotic expansions are used to define a parameter 
regime where a viscous-diffusive-buoyancy balance 
is dominant. The solution, expressed in terms of 
an inverse Fourier transform, is numerically inte- 
grated. The results show that, as in the non-dif- 
fusive problem, a region of closed streamlines 
exists ahead of the body. However, unlike the 
case where diffusion is neglected, the density field 
within this reci rcuiating region is uniquely determined 
and found to be statically stable. It is also found 
that varying the relative amount of diffusion affects 
not only the density distribution, but the velocity 
profile as well, Indicating a strong coupling between 
the vorticity and specie equation. 
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INTRODUCTION 


It Is a well-known feature of stratified fluid motions that the 

influence of a body moving horizontally tan extend far in the upstream 

and downstream directions. In the limit of very large stratification 

and negligible viscous forces, Ylh' has shown that this influence is 

contained in a slug of fluid of infinite horizontal extent. In reality, 

however, viscous and diffusive effects will modify this behavior sc 

2 

that far from the body the fluid is undi * turned. Long has shown that, 
far upstream of an obstacle, the velocity perturbation decays alge- 

- 3 /if 

braically as X as a result of viscous stresses. In a subsequent 
work, Long ^ showed that diffusion of the stratifying agent across 

streamlines alters this far field behavior such that the velocity decays 

- 2/3 - 2/3 4 

as X . This X dependance was also obtained by Koh for the 

case of a viscous-diffusive flow towards a sink. 
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Graebel , using Fourier integral techniques, extended Long's far 

field similarity analysis and found a solution for the upstream influence 
of a body in a nondiffusive fluid which is uniformly valid for X^0(R e R.), 
where R g and R. are the Reynolds number i.nd Richardson number respec- 
tively. Graebel showed that this solution is vclid provided R.>^l, 

R e R i»S and the Prandtl number is infinite. However, Browand and 
Winant^ in discussing Graebel's'* analysis, showed that this solution is 
uniformly valid for the entire flow field, except in the neighborhood of 
the singular points at the edges of the body, provided that the addition- 
al requirement^ R /R.^Xl is satisfied. The calculations of Browand and 
c i 

Winant^ show that immediately ahead of the body there is a large recirculation 
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region containing closed streamlines, and they experimentally verified 
the existence of such a region. However, because diffusion of the 
stratifying agent was neglected, a unique solution for the density dis- 
tribution within this region containing closed streamlines does not exist. 
Thus, Browand and Wlnant^ were unable to present a realistic description 
of the density structure within this region. 

The present study extends Graebei's analysis to include the diffusion 
term In the energy equation. The results show that when the stratifying 
agent is allowed to diffuse across streamlines, the density structure within 
the recirculating region is uniquely determined, and the calculations 
reveal a density structure whirl is physically plausible. The method 
employs asymptotic expansions applied in a straightforward manner, and 
it is shown that the elegant procedure followed by Freund and Meyer^ who 
also considered this problem, is not required. 

ANALYSIS 

The problem to be considered is a linearly stratified infinite fluid 
medium flowing past a two dimension transverse flat plate, as depicted in 
Fig. 1. 

For a steady two d imens iona I flow of an incompressible Boussinesq fluid, 
the equations governing the vorticity and energy are given by: 

L <*• - iar v ] v2 ' 1 ' "> 

L (x, z.f) 
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where, 


1 • njj-r.fc 


g (*•*) - (f -fit ) * f (x,z ) 


r -afiL®. 

1 , ! f 
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. , S ■ ^ , u - f , w ■ - y 
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R - -?■ 


Coordinates are nond i mens iona 1 i zed by h, velocities byU^, and i3andJ are 

respectively the kinematic viscosity and diffusion coefficient. 

The far field asymptotic behavior of Eqs. 1 and 2 may be examined 

by compressing the horizontal coordinate by , Rj. $ c ) such that 

x “ x 2P (/*<!< 1), and scaling the density perturbation by £ (R , R j • S^) 
i — rr i _ 

such that g (x,z) » o J (x,z) ($4<1). By introducing this scaling 
into the governing equations it is found that a viscous-buoyancy balance 
exists in the far field provided 


** " k Ja.s > ** 1 

eV | c 



It is also found that the nonlinear terms in the vorticity and energy 
equation may be neglected provided 



The equations, to this order, may be combined to yield 


, 6 - 


.2 c' 


9z 6 3i- 


( 3 ) 


I'q. 3 may be solved using Fourier transform techniques. By imposing a 
boundedness constraint upon | as x«^ + cu , the solution may be written 
as : 

F(k)e K ~sin(kz)dk (4) 


f ' (T,z) - ^ F 


.3 X 

k sin(kz)dk 


where F (k) Is the Fourier sine transform of the "initial condition" at 
x - 0. The derivation of this quantity is the crucial element in the 
solution of the present problem (and also in the solution presented by 
Graebel"’), but the physical reasoning entering Its determination has not 
here-to-fore been made explicit. Hence, we present a discussion of the 
intricacies involved in both the Graebel'* and Freund and Meyer^ solutions. 

The important elements in specifying the "initial condition" are that 
the horizontal mass flux is conserved and that the integral solution (4), 
when written in terms of Long's^ upstream wake similarity solution (which 
must hold as x — + a> ) , matches uniformly with that solution. That is, the 
correct initial condition must yield a solution satisfying the constraint 
of constant momentum flux defect in *’he upstream wake and continuity of 
horizontal mass flux at x * 0. This is sketched illustratively in Fig. 2. 
The only conditions meeting these two requirements simultaneously are those 
„ lt h ' 

u(x “ 0, z) 

for nondiffusive mot ion' and 


-|$(z + D + £(z - 1)j, 


M (5a) 




u(x - 0,z) ( z 2 _, ) * ^ £j 2 | + ( z 2 - 1 ) ] , | Z |*J 


(5b) 
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for the diffusive case. The first condition was used without comment by 
Graebel"’, and the latter condition was derived via an integral equation 
approach by Fruend and Meyer^. Note that the effect of diffusion is clearly 
seen in spreading the initial condition and weakening the singularities at 
z - +1. The strength of these (integrable) singularities are determined 
solely by the mass r I ux requirement, and their type is determined by the 
asymptotic momentum flux requirement. 

Thus, by imposing upon Eq. b the condition at x = 0 suggested by 
Fruend and Meyer and given by Eq. i>b, we obtain their solution which is 
wr i tten as: 



ao 


(x,z) 



w 


-k 3 ; 


cos(kz) 
k^ sin(kz) 
sin(kz) 


dk ( 6 ) 


cos (kz ) ' 


k sin(kz) 

The solution given by Eq. 6 is valid within the region x = 0 (R^^/rTs^). 

Close to the plate, the equations may be rescaled .hear layer coordinates 

(see Tig. 3 ), Within this region a second solution ; ;ing a viscous-buoyancy 

balance may be obtained by stretching the vertical coordinate as: 

z -1 


(z-1 ) 




1/3 


Nonlinearity may be neglected provided 

1/3 


w 


21 1/3 
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and the resulting equation is identical to Eq. 3, which governs the far 
field flow. However, contrary to the nondiffusive problen (Browand and Wi- 
nant^), no similarity solution Is possible in this case. Thus, under these 
more stringent cond i t ionSj the entire flow field, except in the neighborhood 
of the singular points at z ■ + I, is governed by a linear vi scous- buoyancy 
balance, and the solution given by Eq. 6 is uniformly valid. The impor- 
tance of the inertial terms near the edges of the plate shows clearly the 
need to derive the appropriate "initial condition" via the procedure dis- 
cussed above - a procedure which circumvents solution of the complete 
equations in the neighborhood of these singular points with subsequent 
asymptotic matching to determine F ( k) . 

RESULTS 

The solution given by Eq. 6 was numerically integrated and the results 
are given in Figs. k-S. Fiqure 4 shows the streamline pattern surrounding 
the body, and it is observed that the recirculation region ahead (and 
symmetrically behind) the body exists even wh' i diffusion across stream- 
lines is not negligible. Although not shown in Fig. , it is noted that, 
as in the case considered by Browand and Winant^ an infinite number of 
such recirculating regions exist. It should be renembered, however, that 
the present solution is uniquely determined within these regions, while 
theirs is not. 

Figure 5 shows the horizontal velocity ahead of the body in the plane 

of symmetry. The stagnation point located furthest upstream occurs at 
-1/3 

x ■ -335, and this defines the extent of the recirculating fluid. Also 
shown in this figure is-Long's^ similarity solution which sat isfactoral ly 
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approximates the flow field for x ' ^ 1.4. 

Figure 6 shows the velocity distribution at various upstream locations 
both within and outside of the region of closed streamlines. One can see 
the marked increase in velocity as the singular points located at the 
edges of the place are approached. 

It is interesting to note that the velocity profile contains a single 
inflection point around z ■ 1 and then monotonical ly decays toward zero. 
This is contrasted with the oscillatory decay of the velocity at large 
values of z for the nondiffusive problem 

The density perturbation profiles at various horizontal locations 
are shown in Fig. 7- It should be remembered that the density perturba- 
tion is scaled by the parameter^"s^7ft7” which must be small if the theory 
is to be valid. Thus, the density pioflle through the recirculation region 
is always statically stable, and varies only slightly from the linear 
gradient of the unperturbed fluid. This shows that, for the present model, 
the velocities are so low and diffusion is so Important that, as a fluid 
parcel is convected about, it is never far from being neutrally buoyant. 

The effect of different diffusion rates upon the velocity and density 
profiles is shown in Fig. 8. The two cases shown are $ c ■ .72 and 6.5, 
which roughly correspond to thermal diffusion in air and water respective- 
ly. It is interesting to note that changing the Schmidt number not only 
affects the density distribution, but also has a profound influence upon 
the velocity profile. T his indicates a strong coupling between the vor- 
ticity equation and the energy equation. 
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FIGURE CAPTIONS 


FIG. I. Geometry and coordinate system. 

FIG. 2. Required horizontal velocity at x - 0 for diffusive and 
nond I f fus I ve case . 

FIG. 3* Shear layer coordinate system close to the singular points at 
the edge of the plate. 

FIG. *♦. Stream line pattern showing the recirculating region. 

FIG. 5* Upstream variation of centerline horizontal velocity. 

FIG. 6. Horizontal velocity profile at various uostream locations. 
FIG. 7* Density perturbation profiles at various upsrream locations. 
FIG. 8. Effect of diffusion upon density perturbation and horizontal 
veloci ty prof i les. 
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